Corrections to Scaling are Large for Droplets in Two-Dimensional Spin Glasses 



A. K. Hartmann 

Institut fur Theoretische Physik, Universitdt Gottingen, Bunsenstrafie 9, 37037 Gottingen, Germany 



(N 
O 
O 
(N 
-i— > 
O 

O 

in 

(N 

a 



c3 



i 

a 

o 

o 



> 

oo 
in 
o 

(N 
O 

c3 



a 

o 

o 



X 
S3 



M. A. Moore 

Department of Physics and Astronomy, University of Manchester, Manchester, M13 9PL, United Kingdom 

(Dated: February 1, 2008) 

The energy of a droplet of linear extent I in the droplet theory of spin glasses goes as I for large 
I. It is shown by numerical studies of large droplets in two dimensional systems that this formula 
needs to be modified by the addition of a scaling correction l _UJ in order to accurately describe 
droplet energies at the length scales currently probed in numerical simulations. Using this simple 
modification it is now possible to explain many results which have been found in simulations of three 
dimensional Ising spin glasses with the droplet model. 
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There are two rival theories to describe the ordered 
state of spin glasses: the droplet theory [[[J |J [| and the 
replica symmetry breaking (RSB) theory of Parisi j|, || . 
The droplet picture assumes that the low-temperature 
behavior is governed by droplet-like excitations, where 
excitations of linear spatial extent I typically cost an en- 
ergy of order l e . Thus in the thermodynamic limit the 
excitations which flip a finite fraction of the spins cost an 
infinite amount of energy if 9 > 0. These excitations are 
expected to be compact and their surface has a fractal 
dimension d s < d, where d is the space dimension. Fur- 
thermore it is usually assumed that the energy of differ- 
ent types of excitations, e.g. droplets and domain walls, 
induced by changing the boundary conditions, are de- 
scribed by the same exponent 9. On the other hand the 
RSB theory assumes that there are low-energy droplet 
excitations in which a finite fraction of the spins in the 
system are reversed but which only cost a finite amount of 
energy in the thermodynamic limit. The surface of these 
excitations is space-filling so that the fractal dimension 
of their surface d s = d Q . 

Many numerical studies have been done in an attempt 
to resolve the controversy, but studies in three J?J and 
four dimensional systems |^| are limited in the range of 
systems sizes which can be studied. In two dimensions, 
however, much larger systems can be studied. For this 
dimension, no spin-glass order exists at any finite tem- 
perature ||, but there are still puzzling results concern- 
ing droplet excitation energies about the groundstate. 
The energies of cross-system domain walls pi |ll| and 
droplet-like excitations [jl2|, [l3| are apparently descibed 
by different exponents, which is not compatible with the 
(usual) droplet picture and might be regarded as a ves- 
tige of RSB type behavior. 

In this letter we study numerically droplets in two- 
dimensional spin glasses similar to the droplets studied 
in Ref. [l^, but up to even larger sizes. We show that the 
numerical results mentioned above can be explained by 
a simple scaling correction. For larger droplets we find a 



behavior which is compatible with the droplet prediction, 
i.e. the exponent 9 is the same for droplets and cross- 
system domain-walls. Smaller droplets appear to be de- 
scribed by an effective exponent 9 . We also study how 
the volume of an average droplet behaves with system 
size. We find that when considering just smaller system 
sizes, droplets appear non-compact, but for larger sys- 
tems, crossover to a compact behavior looks likely. Thus 
in two-dimensions the droplet picture seems to be entirely 
valid. Finally, we discuss also some recent numerical re- 
sults [l^, in three dimensions, which also may be 
explained within the droplet picture using the correction 
to scaling for droplet energies introduced here. 

We consider the Hamiltonian which is usually studied 
in numerical work of the Ising spin glass model: 



H 



(1) 



where the sites i lie on the sites of a cubic lattice with 
N = L d sites, Si — ±1, the have a Gaussian distribu- 
tion of zero mean and unit variance and couple nearest- 
neighbor sites on the lattice. 

It is important to realize what the low-energy excita- 
tions look like in spin glasses. Consider a domain wall 
crossing a system M X M X L. The energy cost of this 
domain on the droplet picture will be of order M e and 
the domain wall will be fractal with dimension d s and 
with an area of order M ds . In three or more dimensions 
the interface may have holes through it. Because it is 
fractal the extent of this wandering is of order M. The 
wandering of the interface by an amount of order M af- 
fects the determination of 6>, see Ref. [l?], where it was 
found that the best results were obtained when L 3> M 
as then the interface is not affected by its interactions 
with the ends of the system. Suppose now we have two 
domain walls across the system. Then if their separa- 
tion is large compared to M they will be unaffected by 
the presence of the other. However, if they are closer 
together they interfere with each other and their overall 
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energy will be greater than if the other one were absent. 
This is because domain walls wander to take advantage 
of weak bonds and since the domain walls cannot cross, 
the presence of the second domain wall will preclude the 
first domain wall from cutting through bonds it would 
have cut in the absence of the second wall. As a conse- 
quence of this, domain walls effectively repel each other. 
If they have a separation of order I, then the repulsive en- 
ergy between them would be expected to vary with their 
separation as a power law, . No investigations of u> 
seem to exist in the literature. 

Consider now a droplet of linear extent /. It too will 
have a fractal surface described by d s . Pictures of large 
droplets have now appeared fl5| , and a systematic in- 
vestigation of them is in fTs[] . Their fractal nature en- 
sures that in three or more dimensions that they have 
holes through them, giving them a sponge-like appear- 
ance. However, because the surface of the droplet may 
wander by a distance of order I the energy of the droplet 
E will be modified by its wandering and "collisions" with 
itself to a form which we suppose by analogy with the 
above is 

E = Al B + Br^. (2) 

The term Bl - ^ is a scaling correction to the form of 
the droplet energy at large I and A and B are positive 
constants. We will show that the scaling correction is 
very significant in two dimensions. 

In this dimension 9 is very accurately determined as 
systems of 480 2 can be studied . The exponent associ- 
ated with a single domain wall, that is is approximately 
—0.29 JO, [l?]]. However, studies where, say, 9 is deter- 
mined from the effects of thermally excited droplets (such 
as Monte Carlo simulations of the spin-glass susceptibil- 
ity) yield an apparent 9' close to -0.47 @ This 
discrepancy has long been a puzzle, and has prompted 
suggestions that perhaps different exponents describe do- 
main wall energies and droplet energies |l2|]. However, 
Eq. (^) offers a way of resolving the discrepancy. For 
small values of Z, droplets would apparently have ener- 
gies decreasing as ' when their energies are affected 
by the correction to scaling term, but at large values of 
I the decrease with I will be slower and be as Z - ' 6 ', as 
expected from the conventional droplet approach. 

We have tested this explicitly by studying larger sizes 
than before using an exact ground-state (GS) algorithm 
p9| . For the special case of a planar system without 
magnetic field, e.g. a square lattice with periodic bound- 
ary conditions in at most one direction, there are efficient 
polynomial-time "matching" algorithms pOfl . The basic 
idea is to represent each realization of the disorder by its 
frustrated plaquettes plj . Pairs of frustrated plaquettes 
are connected by paths in the lattice and the weight of 
a path is defined by the sum of the absolute values of 
the coupling constants which are crossed by the path. A 



ground state corresponds to the set of paths with mini- 
mum total weight, such that each frustrated plaquette is 
connected to exactly one other frustrated plaquette. This 
is called a minimum- weight perfect matching. The bonds 
which are crossed by paths connecting the frustrated pla- 
quettes are unsatisfied in the ground state, and all other 
bonds are satisfied. 

We want to mimic the droplet generation used in Ref. 
[l2] , where a center spin was forced to reverse its ori- 
entation with respect to its orientation in the GS, the 
boundary spins were fixed, and a new GS under these 
constraints was calculated. Note that fixing the absolute 
orientation of a spin does not allow for a fast polynomial 
algorithm. Hence, a different approach must be used to 
study large systems. 

Our generation of the droplets works in the follow- 
ing way. For each realization, first a ground state is 
calculated with free boundary conditions in both direc- 
tions. Then several "hard" bonds are introduced, i.e. 
bonds with a high value of the absolute strength (e.g. 
\Jij\ = 2N x max^^jjy}). If the subsystem of hard 
bonds does not exhibit frustration, no hard bond will 
be broken, if a new GS is calculated. First, all bound- 
ary spins are fixed relative to each other by hard bonds 
around the border chosen to be compatible with their 
GS orientations, i.e. the bonds between pairs of bound- 
ary spins are replaced. Second, a line of hard bonds is 
created which runs from the middle of (say) the left bor- 
der to a pre-chosen center spin, again fixing the pair's 
spins in their relative GS orientations. Next, the sign of 
exactly one hard bond on this line is inverted. Finally, 
a GS of the modified realization is calculated. With re- 
spect to the original GS, the result is a minimum energy 
droplet fulfilling the constraints that it contains the cen- 
ter spin, does not run beyond the boundary and that it 
has a surface which runs through the hard bond which 
has been inverted. The energy of the droplet is defined 
as the energy of the resulting configuration calculated us- 
ing the original bond configuration. For each realization, 
this procedure is iterated over all the bonds which are lo- 
cated on the line from the boundary to the center, when 
in each case exactly one hard bond is inverted. Among all 
generated droplets, the one exhibiting the lowest energy 
is selected. 

This procedure is similar, but not equivalent to the 
method used in Ref. The difference is that the surface 
of our droplets crosses the line from the boundary to the 
center spin exactly once. This means, it cannot cross 
this line more than once, i.e. it cannot "meander" across 
this line. But it may meander in all other directions. To 
reduce the influence of this constraint, we have iterated 
over all four lines of bonds running from the left, right, 
top and bottom boundary to the center spin, and selected 
the minimum energy droplet among all the 2L— 2 droplets 
generated in this way. Hence, we have obtained droplets 
which may meander in any direction, but not in all four 
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directions at the same time. 

We have tested explicitly the influence of this con- 
straint by measuring in how many directions each of the 
lowest energy droplets meanders. For small systems me- 
andering does not occur, but meandering increases with 
system size. For the largest systems we have considered, 
about half the droplets do not meander at all, about 33 
% meander in one direction, 14 % meander in two direc- 
tions and three percent in three directions. Hence, for the 
system sizes studied, four-way meandering droplets prob- 
ably would occur rarely, if it were allowed by the method. 
Hence, we believe that our results are not strongly influ- 
enced by the constraint imposed by our algorithm. Fur- 
thermore, we have also tested droplets generated in a 
different way, which allow for full meandering, but this 
approach is restricted to small sizes, see below. 

We have generated droplets for system sizes in the 
range 6 < L < 160, each time 5000 independent realiza- 
tions. The average droplet energy as a function of system 
size is shown in Fig. |l|. For small sizes up to L = 40, the 
energy can be fitted by an algebraic function ~ L e with 
0' = —0.45(1), which is compatible with the earlier result 
jl2| which had been interpreted in favor of the droplet 
exponent being different from the domain- wall exponent. 
For larger sizes, a clear crossover can be observed, the 
data is compatible with an exponent 9 — —0.29. A fit of 
the exponent in [80,160] yields 9 = —0.32(1), which in- 
dicates that the thermodynamic limit may not been fully 
reached. We have also tested the scaling ansatz Eq. (^) 
by fixing 9 = —0.29, resulting inu = 0.97(5) and a very 
good quality p3] Q = 0.81 of the fit. Curiously, the value 
for to is close to the value which 9 takes in one dimension, 
-1. 

To estimate the influence of the slightly restricted me- 
andering, we have also studied droplets, where in ad- 
dition to the hard bonds around the boundary, only 
one additional inverted hard bond was introduced. For 
each system, this inverted hard bond was iterated over 
all bonds connecting at least one non-boundary spin. 
The minimum energy droplet was selected among all 
droplets containing the pre-selected central spin. Due 
to the large numerical effort, only sizes up to L = 30 
could be studied. When considering up to 1000 realiza- 
tions, the data was compatible with an algebraic decay 
with 0' = —0.47(1). We have increased the statistics, 
such that we have treated 80000 realizations for small 
sizes down to 20000 realizations for L = 30. Now the 
straight algebraic fit has a very low quality (Q = 10 -26 ) 
compared to the fit including the correction to scaling 
(tu = -2,9= -0.29, Q = 0.61). 

For further comparison with Ref. [l^, we show in Fig. 
^ the average droplet volume and surface as a function 
of the system size. By fitting to algebraic laws, we find 
a volume fractal dimension of D v = 1.80(1) and a sur- 
face fractal dimension of d s = 1.10(1), which are both 
compatible with the results obtained before. Hence, the 
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FIG. 1: Average droplet energy E as a function of the system 
size. The dashed straight lines represent algebraic functions 
with exponents 9' = -0.47 and 9 = -0.29, while the full 
curved line represents the function from Eq. ([]) with u = 1, 
9 = -0.29). 
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FIG. 2: Average droplet volume V and average droplet sur- 
face S as a function of the system size. The straight lines 
represent algebraic functions with exponents D v = 1.8 re- 
spectively d s = 1.1. 



crossover observed in the energy is not obvious when con- 
sidering volume and surface and the droplets seem to be 
non-compact (i.e. D v < 2). We have tested, whether this 
result may be a consequence of a correction to scaling, 
by fitting V(L) to a function of the form EL 2 + FL D '- 
resulting in D' v = 1.52(5) and a much better quality of 
the fit; Q = 0.59 as opposed to Q — 0.002 for the pure al- 
gebraic fit. Hence a crossover to compact droplets seems 
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likely. In the figure, both functions would be hard to dis- 
tinguish, hence we have omitted including the correction 
to scaling fit. We have also tried a correction to scaling 
for the surface, since d s around 1.3 in other studies ]To| ], 
but using the correction to scaling we could change d s in 
a wide range without affecting the quality of the fit. 

A crossover between large and small L dependencies 
was seen in the work of Middleton Q. He studied the 
scaling of ground-state link overlaps and observed a grad- 
ual evolution of the finite-size behavior. Middleton him- 
self attributed this gradual evolution to averages over 
droplets of all sizes up to the system size L. This kind of 
behavior is an additional mechanism, independent of our 
scaling correction for the droplet energies themselves, to 
explain why certain averages might be slow to converge 
to the expected behavior for large L. 

We have seen that in two dimensions that the correc- 
tions to scaling are large, so that for medium-size sys- 
tems they can mask the true asymptotic behavior. A 
similar effect would be expected to occur in higher di- 
mensions and indeed might already have been observed 
by Lamarcq et al. |l4| in their study of droplets in three 
dimensions. They first computed the ground state for 
both N = 6 3 and N — 10 3 spins. They then chose an 
arbitrary reference spin and flipped it along with a clus- 
ter containing V — 1 other spins connected to it. They 
next minimized the energy of this cluster by exchange 
Monte Carlo, but with the constraint that the reference 
spin is held fixed and the cluster was always connected 
and of size V . They found the largest extension (mean 
end-to-end distance) / of the cluster and found that for 

V < 33 its energy E(l) varied as l e , where the exponent 
9 was —0.13 ± 0.02. For larger V values they found that 
E(l) was increasing rather than decreasing with I but 
they were unsure whether this might not be an artifact 
of insufficient numbers of Monte Carlo steps. Note that 
according to Eq. (||), when both 9 and uj are small and 
positive, E(l) will have a shallow minimum at some value 
lo, similar to the behavior actually seen in Ref. [l4| 

Furthermore, for droplets whose linear dimensions are 
similar to the minimum at lo, the apparent value of 9 
would appear to be close to zero. This may be the origin 
of the behavior in Refs. [l5|, |l6|, where simulations of 
Parisi's spin overlap distribution function and the link 
overlap function were only consistent with the droplet 
picture if 9 was taken to be about zero rather than its 
domain wall value of order 0.20 in three dimensions J^j] 
and 0.70 in four dimensions ||. 

To summarize: there exists a natural correction to the 
usual scaling formula for the energy of a droplet. We 
have explicitly verified its validity by numerical studies 
of large systems in two dimensions. If a similar correc- 
tion is valid in three dimensions, it can explain a number 
of puzzling features which have been seen in simulations. 
A direct check of the validity of this scaling correction 
needs numerical studies of larger three-dimensional sys- 



tems than those currently studied. 
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